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Abstract. The problem of (non)random distribution of points on sphere and in space is 
investigated. The procedure for obtaining preferred direction (and plane) for points on 
the sphere (in the sky) and in the space is discussed. 



1. Introduction 

Klacka (1998) has discussed the methods for finding preferred direction and preferred 
plane used for several decades in application to perihelion concentration of long-period 
comets. The aim of this paper is to present correct methods for finding preferred direction 
and preferred plane for points situated on the sphere, and, also, in general, in the space. 



2. Preferred Direction 



2.1. Points on Sphere 

Points on the sphere may be characterized by direction cosines li, rm, Ui. By definition 
of direction cosines, the relation If + mf + nf = 1 holds. 

If we have N points on the sphere, we have a set of N triples of direction cosines: 
{h}iLi: {"^^}fc=lI {ni}iLi (subscript i denotes the z-th point). 



Let we want to find a preferred direction from the set of N points situated on the 
sphere. Formally, wc can construct a function of the preferred direction characterized by 
the direction cosines I, m, n, or, tIq = ( ^, m, n): 

N 

$K) = 2-i J2w, (no - n,f - A (n^^ - l) , (1) 

where Wi{> 0) is the weight of the i—th point and A represents the multiplicator of La- 
grange. The preferred direction Uq is characterized by the requirement that it minimizes 
the function $ (rio). The result is 

I ^ i'^Wi li\ / fi , m = I ^ Wi m^ I / jj, , n= i'^Wi ni\ / /j, , 



' N \^ / N \^ / N \^ 

\ \'^Wi hj + i'^Winiij + i'^Wi Uij . (2) 

Definition of $ (rio) offers also the measure of spread: MS = ■\/$ (uo) = ^/l — /i, if 
Si=i Wi = 1 (normalization of weights); since < /i < 1, < MS < 1 - the larger the 
value of MS, the worse determination of the preferred direction (the case /i = means 
that Uo docs not exist). Thus, in reality, parameter fi can be used as a measure of the 
significance of the existence of the preferred direction. 

However, the important question is: What is the motivation for the special values of 
the weights Wi and of the form of the function $ {uo)'^ Other forms of functions $ (no) 
would yield other preferred directions tIq. The argument in favour of the choice of $ (no) 
is in physical equation (4) in section 2.2 of this paper (all points arc situated on the 
sphere) . 

2.2. Points in Space 

We have discussed the determination of the preferred direction rio in section 2.1. The 
procedure was defined by the function <i>(no), which determined also the measure of 
spread (we have shown that the quantity /i may be used, equivalently) . This represents 
purely geometrical method. 

Another method is of a physical character. Using any physical quantity / (radius 
vector, velocity vector, etc.), we can minimize the function 



*'(/o)=|e^« ifo - f^A I |e^./.4 



_ (3) 

The result is 

" (4) 



So^l^W^fA I V^Wi\ , 



and, 



$'(/o)=i - |(e-^) fA / ifi^^fA 



(5) 



We have < $' (fg) < 1 for fo given by Eq. (4). Thus, the quantity $' (/o) may serve 
as a measure of spread, or, the measure of the significance of the value fo- In the case of 
kinematical quantities Wi = rrii (masses) and fo represents the corresponding quantity 
for the center of mass. 

The preferred direction is determined by the relation n'^ = fo / \fo\- The better ap- 
proximation of $' (/o) to 1, the less significant is the direction n^. Of course, n^ is 
different from Uo discussed in section 2.1. In the case \fi\ ~ f, i — 1 to N, both meth- 
ods yield the same results (n^ = Uo). Approximately, the unit vector < X]i=i 9i ''^j r / 
1 I X]i=i 9i ''^d f ~ Qi~ perihelion of the i—th comet - corresponds to the unit vector 
< X]i=i ''^i f / 1 I X]i=i '^il f if tl^c distribution in q is independent on the observed di- 
rections n. 



Eq. (5) may be rewritten to the form 



N / ^ \ 



2.3. "Local" Character of the Measure of Spread 



(6) 



Formulae of the measures of spread discussed in sections 2.1 and 2.2 represent "global" 
character. This means that the quantities fi, <!>' will yield a result "direction is relevant" 
already in the case if only one of the components of fo, or, one of the spherical angles 
( a, (5 ), is very significant. Methods of the sections 2.1 and 2.2 cannot determine if all 
components determining preferred direction are significant. 

The problem just discussed may be solved by applying central limit theorem. We 
calculate for this purpose 



joj ± O'foj , O'foj 



N 






(analogously for directional cosines - fox, foy, foz are mean values of directional 
cosines, in this case). Standard tests of mathematical statistics will decide about relevance 
of the quantities fox, foy, foz- Great advantage of this method is that, in general, not 
all of the quantities fox, foy, foz are significant (different from zero). This can help us 
in better understanding of the observational data, in their physical and astronomical 
interpretation. 



3. Preferred Plane 

3.1. Points on Sphere 

3.1.1. Center of the Sphere Lies in the Plane 

Let the preferred plane, eharacterized by the unit normal vector p, is obtained from 
minimization of the function 

N 

e{p)^J2w, {p-Uif - A (p2 _ 1) . (8) 

1=1 

Stationary points are characterized by the equations X]i=i ^« (P ' '^O ''^i ^ ^ P = 
0, O (p) = X = X]i=i '^i iP ' '^i) — 0- There is minimum of O (p) for X{min), maximum 
for X{max) and there is no extreme of O (p) for A corresponding to X{min) < A < X{max) 
(the corresponding quadratic forms are positively definite for A(?nm), negatively definite 
for X{max) and indefinite for the third case ~ only two of the three stationary points 
correspond to extremes). 

Moreover, 

N 

X{min) < (1 /3) ^w^ . (9) 



i=l 



N 



The last condition yields X{min) < 1/3 for the case X]i=i ""^j ~ 1' and, thus, 
X]i=i ""^j (P ' "«) < 1 / 3; in other words, the quantity X]i=i ''^i (P ' "■«) / (1/3) 
compares the observed distribution with respect to the uniform (equal weights) distribu- 
tion. 

The condition X{rain) < (1/3) X]i=i ^« '^^^ ^^ easily proved as follows. O (p) = 
'^i^i'Wi (p ■ rii) we write as & (p) = p^ F p, where the tensor F is defined by its 
components Tjk = X]i=i ^i {''^i)j {f^i)k- '^^'^ symmetry of F enables us to find an or- 
thonormal base in which F is diagonal (Tr(F) is invariant of the transformation orthog- 
onal transformation; primed components), i. e., Tjk = X]i=i "^i (''^0 ^jk- The equation 
r p = Xp yields for the eigenvalues Ai = X]i=i ^i ('^Oi' ^^2 = J2t^i wi {n[)^, A3 = 
Y.i=i 'Wi ("•O3' i- *5-' Ai + A2 + A3 = J2i=i '^i^ and, thus, X{min) < ( 1 / 3 ) Y.i=i ^*- 
The relation Ai + A2 + A3 = X]i=i ^i niay serve as a control in practical calculations. 

3.1.2. General Plane and Sphere 

If the preferred plane is obtained from minimization of the function 

N 

Q{p,d)^Y.'^^ (P-^* + '^f - A (p2 _ 1) , (10) 



then the process 

of minimization yields for the minimum: X]i=i ^« {p ' f^i + d) rii — X{min) p = 0, 
J2i=i Wi {p-rii + d) = 0, Q (p, d) = X(min) = J2i=i ^t {P ■ ^i + df > 0. 

Moreover, \{min) < (1/3) S X]i=i ^i ~ (X]i=iWi) M^ f; where ^ = 



2 



(X]i=i^i^i) + (Si=i^j"^i) + ( X]i=i '^« '^i ) (tttc same quantity as in sec- 
tion 2). The last condition yields X(min) < (1 / 3) (1 — /i^) for the case X]i=i ^i — 
1, and, thus, X]i=i ^i {p ' f^i) < (1 / 3) (1 — /i^); in other words, the quantity 
X]i=i ^i (p • rii) / (1/3) compares the observed distribution with respect to the 
uniform (equal weights; ii — 0) distribution - Eq. (9) may be used. 

The condition X{min) < (1 / 3) < X]i=i ^i ^ ( X]i=i ^0 M^ r *^^^ t)e easily proved 
as follows. 6(p, rf) = X]i=i ^i {p ■ rii + d) — A (p^ — l), we write as Q{p,d) = 
p^ Fp + 2 d p^ a + d^ (X]i=i^i) ^ '^ {p^ P ^ l)j where the tensor F is defined 



the eigenvalues A^ = ^^^^ w, (nQ^^ - (x;«^i ^*) E^^i ^* «)/ 

-^ ,„.^ f\^^ ..,.\^^ ,,2 



by its components Tjk = X]i=i '^i ('^Oi (''^Ofci '^.nd the vector a = X]i=i (^i '^O- 
The conditions of stationarity yield Fp + da = X p, p^ a + d ( X]i=i wA = Q . 
The last two equations yield F p — I X]i=i ^0 (p • a) a = A p, or, in components 

^jk ^ \l2ii=i'^i) ^j o,k(Pk = ^Pjy or, shortly Hjk Pk ~ ^ Pj (summation over 
repeated indices is supposed in the last two equations). The symmetry of Ft enables us to 
find an orthonormal base in which Ft is diagonal (Tr(n) is invariant of the transformation 
- orthogonal transformation; primed components). The equation Hp = Ap yields for 

J = 1, 2, 3, i. 
e., Ai + A2 + A3 = [J2'i=i'^i) ~ (Z^iLi'^i) M^> where fi = \a\. Thus, X{min) < 
(1 / 3) ^ X]i=i ^i ^ ( X]i=i ^0 M^ f- The relation for the sum of A— s may serve as a 
control in practical calculations. 

Equation X]i=i ^i {p ■ rii + d) = and the results of the section 2 yield |d| < 1. 

Example : Fet us consider {n,}^^^ = { (1, 0, 1) / \/2; (0, 1, 1) / \/2; (- 1, 0, 1) 
/ V2; (0, - 1, 1) / V2 } ; w, = 1/4, « = 1 to 4. 

The method described in section 2 yields for the preferred direction: Uo = (0, 0, 1), 
fx = V2/2 = 0.707. 

As for the preferred planes: 
i) e (p) = X{min) = 1 / 4, p = (pi, p2, 0), pj + p^ = 1. 
u) e (p, d) = X{min) = 0, p = (0, 0, (±) l),d = - (±) V2/2. 

Particular results: a = (0,0, l/\/2), X{min) < 1/6, F = dia5(l/4, 1/4, 1/2), H = 
dm5(l/4, 1/4,0). 



3.2. Points in Space 

3.2.1. Center of the Coordinate Axes Lies in the Plane 

Let the preferred plane is obtained from minimization of the function 

N 

G{p) = Y,^^ {p-f^f - A (p2 _ 1) , (11) 

where fi is "radius" vector of the z— th point in the (phase-)space. Stationary points 
are characterized by the equations X]i=i ^i {P ' fi) fi ^ X p ~ 0, O (p) — \ — 
X]i=i ^i ip ' fi) ^ 0- There is minimum of Q {p) for X(min), maximum for X{max) 
and there is no extreme of O (p) for A corresponding to X{min) < X < X{max) (the 
corresponding quadratic forms are positively definite for X{min), negatively definite for 
X{max) and indefinite for the third case - only two of the three stationary points corre- 
spond to extremes). 

Moreover, 

1 ^ 
A (mm) < - ^ w» /^^ . (12) 

i=l 

The last condition yields X{min) < X]i=i /i^/(3^) for the case Wi — l/N, and, thus, 
Eili Wi (p • /j)^ < Eili fi'^/i^N); in other words, the quantity 3 Y.f=i (P ' Z^)^ / 
< X]i=i fi^ ( compares the observed distribution with respect to the uniform distribution. 

The condition X{min) < (1/3) X]i=i ^i fi^ c^'^ t)e easily proved in an analogous 
way as it is presented in section 3.1.1. 

3.2.2. General Plane and Space 

If the preferred plane is obtained from minimization of the function 

TV 

e(p,d) = ^u;, (p-/, + df - A (p2 _ 1) , (13) 

i=l 

then the process of minimization yields for the minimum: X]j=i ^i {P fi + d) fi — 
X{min) p = 0, Yln=i ^i (P ■ fi + d) = 0, e (p, d) = X{min) = J2iLi '^i ip ■ fi + df . 

Moreover, X{min) < {J2iLi "^^ f%^ ~ yl^=\ ^i) a^}/3, where a = X^ili (^^ /O- 
The last condition yields X(min) < { [ (l/N) X]i=i fi^ ] ^ ^^ }/3 for the case Wi = l/N, 
and, thus, J2i=i Wi {p ■ fi + d) < {[ {l/N) Y,i=i fi' \ ^ °^ }/3; in other words, the 
quantity X]i=i {P ' fi + d) / { (1 / 3) J2i=i fi^ } compares the observed distribution 
with respect to the uniform (a = 0) distribution - Eq. 12 may be used. 

The condition X(rain) < (1/3) {J2i=i ^i fi^ ^ (Si^i^i) ^^} can be easily 
proved in a manner analogous to that presented in section 3.1.2. 



4. Simultaneous Consideration of Preferred Direction and Preferred Plane 

We present a general discussion for simultaneous calculation of the preferred direction 
and the preferred plane. 

Let us consider a function 

N N 

$K,p) ^ 2-1 ^u;, {no - n,f + A^w, {p ■ mf - A {uo^ - l) (14) 

1/ (p-rio) , 

which is minimized in order to obtain the preferred direction rio and the unit vector p 
normal to the preferred plane containing the center of the unit sphere; Wi > 0. Coefficients 
A, p, V are multiplicators of Lagrange. The coefficient A (> 0) is arbitrary and there is 
no argument which value of A one should take into account. The case A = corresponds 
to the case $ (rio) discussed already in section 2. The case X]i=i ^i ('"■o — ni) <C 
A X]i=i ^i ip ' f^i) ^-iid A = i^ = corresponds to the case when there is no apriori 
relation between tIq and p. 

We present two simple examples. 
1- {r^^}tl = {(1-0, 0); (-1, 0, 0); (0, 0, 1)} ; w, ^ 1/3, * = 1 to 3. 
It can be easily verified that Uo = (0, 0, 1), p = (0, ±1, 0). Thus, Uo ■ p = 0- 
Moreover, in order of presenting some other properties, we offer some other results: 
< X]i=i '^i (■'^o • rii) \ / \ X]i=i 'Wi\ = 1/3; /i (defined in section 2) = 1/3, probability 
that jj, > 0.33 in the case of random uniform distribution is 83 % ; -^ X]i=i ^i (P ' '^O r 
/ < X]i=i Wi> = = \{min); definition wi = W2 = l/(fc + 2), W3 — k/(k + 2) yields 
\ J2i=i ^i ('^o • ni) > I \ X]i=i ^i f = kjik + 2) and this may be greater than 1/3 for 
fc > 1, ^(fc) = kjik + 2), X(min, k) = 0. 

2.K}ti -{(l,0,0);(-l,0,0);(0,l,0);(0,-l,0);(0,0,l)};w;, = 1/5, z-lto5. 
It can be easily verified that n,, = (0, 0, 1), p = (0, 0, ±1). Thus, Uo ■ p = ±1. 

The last example shows that the function ^{no,p) is reasonable only in the case 
v = 0. Thus, the important property of Eq. (14) is that minimization of $ {no,p) yields 
independent equations for Uq, p. 

5. Physics and Geometry 

Finally, we show the relation between physics and geometry. We have 0' (p) = 
X]i=i ^i (P ' fi) — 12i=i'^i l/iP (P ■ ■'^i) , which is equivalent to O (p) for the case 
Wt = w'^ l/ip, i = 1 to iV. The minimization of $' (/o) = YlfLi'^i ifo - ftf 
yields /o = {Eti< /.} / {Y^=i<] - {Eti^ m "4 / {T^^iK] , "o = 



{ YliLi w'ilfil m\ / 1 1 Y^iLi w[ \fi\ nj| \ . This is equivalent to the minimum of $ (Uo) 
= Hi=i'^i ("o - riif for the case wi = w[ |/i|, i = 1 to A^ (or, wi = w- |/j| / 

{lEti^l/dnd})- 
6. Conclusion 

We have presented several methods for obtaining preferred direction and preferred plane. 
Methods determining significance of the obtained results are also presented, and, also, 
tests on uniform distribution can be easily elaborated using computer modelling. 
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